Parity- Affected Persistent Currents in Superconducting Nanorings 
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We argue that at sufficiently low temperatures T superconducting parity effect may strongly 
influence equilibrium persistent currents (PC) in isolated superconducting nanorings containing 
a weak link with few conducting modes. An odd electron, being added to the ring, occupies the 
lowest available Andreev state and produces a countercurrent circulating inside the ring. For a single 
channel quantum point contact at T = this countercurrent exactly compensates the supercurrent 
I e produced by all other electrons and, hence, yields complete blocking of PC for any value of the 
external magnetic flux. In superconducting nanorings with embedded normal metal the odd electron 
countercurrent can "overcompensate" I e and a novel "7r/7V-junction" state occurs in the system. 
Changing the electron parity number from even to odd results in spontaneous supercurrent in the 
ground state of such rings without any externally applied magnetic flux. 

PACS numbers: 74.78.Na, 73.23.Ra, 74.45. +c, 74.50.+r 
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Thermodynamic properties of isolated superconduct- 
ing systems are sensitive to the parity of the total num- 
ber of electrons Q, 0] even though this number J\f is 
macroscopically large. This parity effect is a fundamen- 
tal property of a superconducting ground state described 
by the condensate of Cooper pairs. The number of elec- 
trons in the condensate is necessarily even, hence, for 
odd jV at least one electron remains unpaired having an 
extra energy equal to the superconducting energy gap 
A. This effect makes thermodynamic properties of the 
ground states with even and odd Af differ. Clear evidence 
for such parity effect was demonstrated experimentally in 
small superconducting islands 0, Q . 

Can the supercurrent be affected by this parity ef- 
fect? In many structures the answer to this question 
is negative because of the fundamental uncertainty rela- 
tion SAfSip > 1. Should the electron number M be fixed, 
fluctuations of the superconducting phase cp become large 
disrupting the supercurrent. On the other hand, in trans- 
port experiments with fluctuations of ip being suppressed 
the parity effect cannot be observed because of large fluc- 
tuations of Af. 

Nonetheless, the parity effect can coexist with non- 
vanishing supercurrent. The way out is to consider sys- 
tems supporting circular persistent currents (PC), such 
as, e.g., isolated superconducting rings pierced by the 
magnetic flux <E>. In accordance with the number-phase 
uncertainty relation the global superconducting phase of 
the ring fluctuates strongly in this case, however these 
fluctuations are decoupled from the supercurrent and 
therefore cannot influence the latter. The aim of this 
paper is to demonstrate that the parity effect may sub- 
stantially modify PC in superconducting nanorings, in 
particular if the total number of electrons is odd. 

Consider, for instance, a single mode quantum point 
contact (QPC) embedded in a superconducting ring. 



Provided the number of electrons in the system is not 
fixed the Joscphson current across this contact is given 
by the well known expression jj] 



m = 



2e ds(ip) tanh eM 



(1) 



ti dtp 2T 
where ip is the phase difference across QPC, e(ip) = 



Ay 1 — Tsm 2 (ip/2) and T is the contact transmission. 
This result has a very transparent physical interpreta- 
tion being related to contributions of discrete Andreev 
energy states E± = ±e(<ys) inside QPC via 



2e 

Ti 



(2) 



dip dtp 

Using the Fermi filling factors for these states f±(E±) = 
[1 + exp(±e(tp) /T)}^ 1 one arrives at Eq. (JTJ. 

Let us now fix the number of electrons inside the ring 
and consider the limit T — > 0. For even Af all electrons 
are paired occupying available states with energies be- 
low the Fermi level (see Fig. la). In this case one has 
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FIG. 1: Andreev levels inside QPC and their occupation at 
T = for even (a) and odd (b) ensembles. 
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f-(E-) — 1, f + (E + ) — 0, the current is entirely deter- 
mined by the contribution of the quasiparticle state E_ 
and Eq. (J2J yields the same result as one for the grand 
canonical ensemble. By contrast, in the case of odd Af 
one electron always remains unpaired and occupies the 
lowest available energy state - in our case E + - above the 
Fermi level. Hence, for odd Af one has f±(E±) = 1 (Fig. 
lb), the contributions of the two Andreev states in Eq. 
(J2J exactly cancel each other, and the current across QPC 
remains zero for any tp or the magnetic flux $. Thus, we 
predict a novel mesoscopic effect - parity-induced block- 
ing of PC in superconducting nanorings with embedded 
QPC @. 

It is important to emphasize that the above physical 
picture is based on the assumption that the supercur- 
rent is determined only by the contributions of two dis- 
crete Andreev levels E±(tp). This is the case provided (a) 
quasiparticle states of the continuum do not contribute 
to the supercurrent and (b) there are no more discrete 
Andreev levels inside the junction. Both these conditions 
are met only for symmetric and extremely short QPC. In 
a general case of SNS junctions with a non-vanishing 
thickness of the normal layer d =/= these conditions are 
violated, and no exact compensation of the supercurrent 
by the odd electron countercurrent is anymore possible. 

Further analysis reveals even richer physics in the lat- 
ter case. In particular, below we will demonstrate that 
at low T and for odd number of electrons a novel ir/N- 
junction state should occur in such SNS rings. This ef- 
fect in turn leads to spontaneous supercurrent flowing in 
the ring in its ground state. Under certain conditions the 
magnitude of this current can be as high as I sp ~ eA/h 
reaching the values up to 10 4- 100 nA for generic BCS su- 
perconductors. Spontaneous currents of this magnitude 
can reliably be detected in modern experiments. 

Parity projection formalism. In order to systematically 
investigate an interplay between the parity effect and PC 
in superconducting nanorings we will employ the parity 
projection formalism 0i IS @ • 

The grand canonical partition function Z(T,fi) — 
Tre -i3(H-^) ig ii n k ec i t0 tne canonical one Z(T,Af) as 



(3) 



Here and below H, is the system Hamiltonian and (3 
I 7 . Inverting this relation and defining the canonical 
partition functions Z e and Z a respectively for even (Af = 
Af e ) and odd (Af = Af ) ensembles, one gets 



Z e/o (T) = 



2vr 



due 



l Z e/o (T,iTu), (4) 



where 



Z e/o (T, n) = i T r { [1 ± (-1)"] e-W-xW } 



-(Z(T,ri±Z(T,» + inT)) 



(5) 



are the parity projected grand canonical partition func- 
tions. For Af 3> 1 it is sufficient to evaluate the integral 
in J3J within the saddle point approximation 



Z e/o(T) ~ e 



-f3(n c/o -^ /o M e/o ) 



(6) 



where f2 e /o — — Th\Z e / (T, n) are the parity projected 
thermodynamic potentials, 



^e/o = % - T In 



i^l± e -/3(^-%)^ 



(7) 



and Q f/b = -Tin [TrKil)^-^^-^)}]. "Chemical 
potentials" n e i are defined by the saddle point condition 
■A4/o = -dn e/o (T,fi e/o )/dfi e/o . 

The main advantage of the above formalism is that it 
allows to express the canonical partition functions and 
thermodynamical potentials in terms of the parity pro- 
jected grand canonical ones thereby enormously simpli- 
fying the whole calculation. We further note that Qf is 
just the standard grand canonical thermodynamic poten- 
tial and fib represents the corresponding potential linked 
to the partition function Z(T, /i + iirT). It is easy to see 
[S| that in order to recover this function one can evalu- 
ate the true grand canonical partition function Z(T, fi), 
express the result as a sum over the Fermi Matsubara 
frequencies ujf — 2irT(l + 1/2) and then substitute the 
Bose Matsubara frequencies ojb — 2nTl (I = 0, ±1,...) 
instead of loj. This procedure automatically yields the 
correct expression for Z(T,fi + iirT) and, hence, for f2/,. 

Having found the thermodynamic potentials for the 
even and odd ensembles one can easily determine the 
equilibrium current I e / a . Consider, as before, isolated 
superconducting rings pierced by the magnetic flux $. 
Making use of the above expressions one finds PC circu- 
lating inside the ring: 



l e/o 



T ± h ~ If 

f e 0(n b -n f ) ±1 - 



(8) 



where the upper/lower sign corresponds to the even/odd 
ensemble and we have defined 



'e/o 



on 



e/o 



9$ 



l f/b 



Homogeneous superconducting rings. Turning to concrete 
calculations we first treat homogeneous nanorings with 
cross section s and perimeter L — 2nR. Throughout this 
paper rings will be assumed sufficiently thin, ^/s < Xl, 
where Xl is the London penetration length. Supercon- 
ducting properties of such rings will be described within 
the (parity projected) mean field BCS theory. At low 
temperatures this description is justified provided quan- 
tum fluctuations of the order parameter, quantum phase 
slips (QPS) 01 1 can be neglected. This requirement can 



be fulfilled only provided the total number of conducting 
channels in the ring remains large ^lj -Wr ~ Pf s ^ 1- 
In addition, the perimeter L should not be too lar ge in 
order to disregard the QPS-induced reduction of PC 
Finally, we will neglect the difference between the mean 
field values of the BCS order parameter for the even and 
odd ensembles 0, • This is legitimate provided the ring 
volume is large enough, V = Ls ^> 1/VA, where v is the 
density of states at the Fermi level and A is the BCS 
order parameter for a bulk superconductor at T = 0. 

The task at hand is now to evaluate the thermody- 
namic potentials fi//&. Within the mean field approach 
these quantities can be expressed in terms of the exci- 
tation energies £k and the order parameter A(r). One 
finds 



Qf = Q - 27 ^ In (2 cosh — ^ 



27V 



Q b = n - 2T^ln(2sinh||) 



(9) 
(10) 



where Q = J d 3 r\ A(r)\ 2 /g + Tr{£}, g is the BCS coupling 
constant and £ is the single-particle energy operator: 

i=-L(-m^---A(r)) +U(r)-fi, (II) 
2m \ or c J 

A(r) is the vector potential and U(r) describes the po- 
tential profile due to disorder and interfaces. 
The excitation spectrum Ek has the form 



£k = e(p) = pv s + \/£ 2 + A 2 , 



(12) 



where p is a quasiparticle momentum, £ = (p 2 — fx) /2m, 
jU = p($) — mv s 2 /2 and 



va = mm„ I n 

2mR V $0 



(13) 



is the superconducting velocity. Both e(p) and vs are 
periodic functions of the flux $ with the period equal to 
the superconducting flux quantum $0 = hc/2e. 

Consider the most interesting case T <C hvp/L (l4|. 
Making use of the above expressions one obtains 

I e = ev S Q e s, I = ev S Q s - e— sgnu s , (14) 

where g e / a = Af e / /V. The first Eq. {H} - together 
with 1)13(1 - coincides with that for the grand canon- 
ical ensemble. In contrast, for odd ensembles (second 
Eq. (|14p) there exists an additional flux-dependent term 
which, however, remains rather small in multichannel 
rings |l5j . Estimating the main contribution to I e / as 
I ~ ev F N r /L, one gets (7 e - I )/I ~ 1/N T < 1. 

Now we proceed to the situations in which parity ef- 
fect gains more importance and leads to new physical 
effects. Namely, we consider superconducting rings with 





FIG. 2: Superconducting ring with embedded SNS junction 
of length d. 



N t S> 1 interrupted by a weak link with only few con- 
ducting channels N. In such systems the mean field BCS 
description remains applicable and, on the other hand, 
the parity effect can be large due to the condition N ~ 1. 
Quantum point contacts (QPC) or, more generally, SNS 
junctions (of length d) can be used for practical realiza- 
tion of such weak links. The corresponding structure is 
depicted in Fig. 2. 

Superconducting rings with QPC. As a first example we 
consider a ring interrupted by QPC. The optimal value 
of the phase difference <p> across QPC is fixed by the equa- 
tion d£l e / /dtp = 0, which yields the current conservation 
condition 7 (c ) = I (r) . Here 7 (c) ~ 2eT A/h is the current 
through the contact and ~ (ev F N r / L)(ip — 27r$/$ ) 
is the supercurrent inside the ring. From the current 
conservation one easily finds [ljj 



ip ~ 27r$/$ , if 7<i*, 
p ~ 27m, if L > L*, 



(15) 
(16) 



where L* = ^qN t /T 3> £0 ~ Kvf/A. In a more general 
case of QPC with N conducting channels in the expres- 
sion for L* one should setT^^f 7^. 

In the limit (|16[1 one recovers the results identical to 
ones for homogeneous rings. Below we will concentrate 
on the opposite limit L <C L* and also assume N <C N z . 
Due to Eq. I|15[) in this case the dependence / e / ($) is 
fully determined by the current-phase relation for QPC 
which can be found by means of Eq. JHJl. Expressing I f n. 
via ip it is convenient to employ the general formula |l7j 



2e 



N 



l f/b 



sm tp 



i=i 



cos^ + Wi(u)f/ b ) 



(17) 



For QPC one has W t {w) = (2/70(1 + h 2 u 2 /A 2 ) - 1. 
Substituting this function into Q17[) and summing over 
Uf one rederives the standard result [4|, while the same 
summation over u>b yields 



.^y^Mcoth^M. 

Ti ^ dip 27 

1=1 



(18) 



where e t {p) = Ay 1 - % sin 2 (<^/2). Finally, the differ- 
ence Qb — flf = flbf is evaluated as a sum of the ring 
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(Q^) and QPC contribut ions. The latter is found 

by integrating If /^(tp) over the phase tp with the result 

N 



Qly = 2T^lncoth 



2T 



(19) 



while the former is defined by the standard expression 
2V / ** bx(ooth £ M) ^.vVATe"^. 



(2tt?i) 3 "V 2T 
Combining all these results with Eq. © we get 



N 



2e^ 



2T 



1 ± 



(coth^) 2 - 1 



2T 



e^/ ]l(coth 

3=1 



A* 



fiMtf -u 1 
2T / 



(20) 



The term in the square brackets accounts for the par- 
ity effect in our system. For N = 1 and at T — this 
term reduces to unity for even ensembles and to zero for 
odd ones, thus confirming our intuitive picture of parity- 
induced blocking of PC discussed above. For T > Eq. 
(|20|l demonstrates that both for even and especially for 
odd JS[ the current-phase relation for QPC may substan- 
tially deviate from that derived for the grand canonical 
ensemble Q, see Fig. 3. For even ensembles the super- 
current increases above its grand canonical value. This 
effect is mainly pronounced for phases tp not very far 
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FIG. 3: The ratio between canonical and grand canonical val- 
ues of PC I e / /Ig.c. (represented by the term in the square 
brackets in Eq. i'201 ) versus tp in a single mode QPC at differ- 
ent temperatures for even (three upper curves) and odd (three 
lower curves) ensembles. Here we have chosen Tg = 0.1A and 
T = 0.99. 



from tp = 7T and - at sufficiently low T - it becomes pro- 
gressively more important with increasing temperature. 
On the contrary, for odd ensembles the supercurrent is 
always suppressed below its grand canonical value. This 
suppression is gradually lifted with increasing temper- 
ature, though at phases tp in the vicinity of the point 
tp = 7r blocking of PC may persist up to sufficiently high 
T. Eq. H2()J) also shows that in QPC with several con- 
ducting channels and at T — > the current through the 
most transparent channel will be blocked by the odd elec- 
tron. Hence, though blocking of PC remains incomplete 
in this case, it may nevertheless be important also for 
QPC with N > 1. 

SNS rings. Finally, we turn to superconducting rings 
containing a piece of a normal metal with non-zero thick- 
ness d. In contrast to the situation of QPC considered 
above, the Josephson current in SNS structures can- 
not anymore be attributed only to the discrete Andreev 
states inside a weak link, and an additional contribution 
from the states in the continuum should also be taken 
into account. Furthermore, for any non-zero d there are 
always more than two discrete Andreev levels in the sys- 
tem. Accordingly, significant modifications in the phys- 
ical picture of the parity effect in SNS rings can be ex- 
pected. 

The key difference can be understood already by com- 
paring the typical structure of discrete Andreev levels in 
SNS junctions (Fig. 4) with that of QPC (Fig. 1). As 
before, in the limit T — > all states below (above) the 
Fermi level are occupied (empty) provided the total num- 
ber of electrons in the system is even (Fig. 4a). If, on the 
other hand, this number is odd the lowest Andreev state 
above the Fermi energy is occupied as well (Fig. 4b) thus 
providing an additional contribution to the Josephson 
current. This contribution, however, cancels only that of 
a symmetric Andreev level below the Fermi energy, while 
the contributions of all other occupied Andreev levels and 
of the continuum states remain uncompensated. Hence, 
unlike in the QPC limit displayed in Fig. 1, in SNS rings 




FIG. 4: Andreev levels in a single mode SNS junction with 
d = 6?ivf/A and their occupation at T — for even (a) and 
odd (b) ensembles. 
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one should not anymore expect the effect of PC blocking 
by the odd electron, but rather some other non-trivial 
features of the parity effect. 

In order to construct a quantitative description of the 
parity effect in SNS rings we will again make use of 
the parity projection technique formulated above. Here 
we only restrict our attention to transparent SNS junc- 
tions where the effects under consideration are most pro- 
nounced. In this case the function Wi(ui) = W(ui) is the 
same for all channels. It reads 



W(w) 



2h 2 



1 ] cosh 



(2ujd 




flLO 2 , [2u>d\ 



(21) 



Substituting this function into (|17J) and repeating the 
whole calculation as above, we arrive at the final result 
which cannot in general be represented in a tractable 
analytic form. Significant simplifications occur in the 
most interesting limit [la T — > 0, in which case we obtain 



eAN ( . tp 

sin — 

h \ 2 



2y sin tp 1 
In — 



T T eA ( ■ ^ 

I = I e - — {sm- 



ysgnip cos tp 



for short SNS junctions y = dA/Tivp <C 1 and 

evpN 



evpN 

Ie = -J-f, 

ird 



/ _ 7rsgny \ 
V N ) 



(22) 
(23) 

(24) 



for long ones d 3> £o ~ Tlvf/A. These results apply for 
— 7r < tp < tt and should be 27r-periodically continued 
otherwise. The term containing ln(l/y) in Eq. (|23|l for 
I e is written with the logarithmic accuracy and is valid 
for ip not too close to ip = ±7r. 



Let us concentrate on the results (|24fl applicable for 
sufficiently long SNS junctions. We observe that at 
T = the current I e - as in the QPC case - coincides 
with that for the grand canonical ensembles Hi, while 
for odd Af the current-phase relation is shifted by the 
value tt/N. This shift has a simple interpretation being 
related to the odd electron contribution (2e/K)dEi/ dtp 
from the lowest (above the Fermi level) Andreev state 
E\(tp) inside the SNS junction. As we have expected, 
this contribution indeed does not compensate for the cur- 
rent from other quasiparticle states. Rather it provides 
a possibility for a parity-induced 7r-junction state |20| in 
our system: According to Eq. I|24l) for single mode SNS 
junctions the "saw tooth" current-phase relation will be 
shifted exactly by tt, see Fig. 5. More generally, we can 
talk about a novel tt/N-) unction state, because in the 
odd case the minimum Josephson energy (zero current) 
state is reached at tp = ±tt/N, see Fig. 6. For any 
N > 1 this is a twofold degenerate state within the in- 
terval — tt < ip < tt 21]. In the particular case N = 2 the 
current-phase relation I {ip) turns 7r-periodic. 

Let us recall that the 7r-junction state can be real- 
ized in SNS structures by driving the electron distri- 
bution function in the contact area out of equilibrium 
H3 H3 Here, in contrast, the situation of a tt- or 
tt/N-) unction is achieved in thermodynamic equilibrium. 
Despite this drastic difference, there also exists a certain 
physical similarity between the effects discussed here and 
in Refs. |22j, |23j, |24{ In both cases the electron distribu- 
tion function in the weak link deviates substantially from 
the Fermi function. It is this deviation which is respon- 
sible for the appearance of the 7r-junction state in both 
physical situations. 

It is also important to emphasize that the countercur- 
rent produced by the odd electron causes a jump on the 
current-phase dependence at ip — 0, and the direction 
of I Q at sufficiently small tp is always opposite to that 
of I e . This is a generic feature of odd ensembles which 



e/o 



N=1 




FIG. 5: The zero temperature current-phase dependence (12 H 
for SNS rings with N = 1: I e (tp) (dashed) line and I (<p) 
(solid) line. 







'odd 










N=4 












N=3 












N=2// 




//■In 


// / -TT 


/// 




//It, 














— ^ 



FIG. 6: The same as in Fig. 5 only for the odd ensembles 
(second Eq. pty l and for N = 2, 3 and 4. 
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FIG. 7: The zero temperature current-phase relation Io{f) 
(—tv < if < 7r) for N = 1 and different values of the parameter 
y = dA/hvF- 



persists for all non-zero values of y. The current-phase 
relation I (y>) for arbitrary values of the parameter y can 
be computed numerically. Examples are presented in Fig. 
7. One clearly observes the current jump at <p = 0. As it 
is obvious from Eq. I|23[l . this feature disappears only in 
the QPC limit y -> 0. 

BCS ground state with spontaneous current. Perhaps 
the most spectacular physical consequence of this cur- 
rent jump is the presence of spontaneous supercurrents 
in the ground state of SNS rings with odd number of 
electrons. Similarly to the case of standard 7r-junctions 
|2flj | such spontaneous supercurrents should flow even in 
the absence of an externally applied magnetic flux. Un- 
like in Ref. |2(], however, here the spontaneous current 
state occurs for any inductance of the ring because of 
the non-sinusoidal dependence / (</?)• 

Consider, for instance, the limit d ^> £o- In the case 
of odd number of electrons the ground state energy of an 
SNS ring can be written in a simple form 



E = 



$2 

2c£ 



$ - 



$ sgn$ 
27V 



(25) 



where the first term is the magnetic energy of the ring (£ 
is the ring inductance and c is the speed of light) while 
the second term represents the Josephson energy of the 
SNS junction. Minimizing Ij25(l with respect to the flux 
$ one immediately concludes that the ground state of the 
ring is a twofold degenerate state with a non-vanishing 
spontaneous current 



evF 



1 



2ev F N £ 



(26) 



flowing either clockwise or counterclockwise. In the limit 
of small inductances C — » this current does not vanish 
and its amplitude just reduces to that of the odd electron 
current at tp — > 0. Hence, at small ring inductances the 



FIG. 8: The spontaneous current amplitude I sp as a function 
of the parameter y at T = 0. In the inset, the same function 
is shown on the log — log scale. Dashed lines indicate the 
asymptotic behavior of I sp (y) in the limits of small and large 

y- 



magnitude I sp of this spontaneous current can easily be 
obtained from Eqs. (|23l24fl : 



= eA 2 d/h 2 VF, 
I sp = evp I ird, 



if d<£o, 
if d>Co- 



(27) 
(28) 



For intermediate values of the parameter y the amplitude 
of the current I sp can be evaluated numerically. The 
results are displayed in Fig. 8. One observes that - in 
agreement with Eq. (|27[) - I sp increases linearly with d 
at small d, reaches its maximum value / max ~ 0AeA/h 
at d ~ £o and then decreases with further increase of 
d approaching the dependence in the limit of large 
d. For generic BCS superconductors the magnitude of 
this maximum current can be estimated as / max ~ 10 -j- 
100 nA. These values might be considered as surprisingly 
large ones having in mind that this current is associated 
with only one Andreev electron state. 

In summary, new physical effects emerge from an inter- 
play between the electron parity number and persistent 
currents in superconducting nanorings. Perhaps the most 
striking observation is that the BCS ground state of such 
rings with the odd number of electrons is the state with 
non-zero spontaneous supercurrent. This and other novel 
features predicted here can be directly tested in modern 
experiments. 

This work is part of the Kompetenznetz "Funktioncllc 
Nanostructuren" supported by the Landestiftung Baden- 
Wurttemberg gGmbH and of the STReP project "Ultra- 
ID" supported by the EU. 
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